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ABSTRACT 


irom awciSe IS uecOncemedsaclth the use ot radio-fre- 
@ulenmey (rt) Current in producing llesions in the brain. The 
moenmeduc Consists Of passing the current between an active 
preetrode implanted in the brain anc an indifferent elec- 
i[mode located on the skull. The resistivity of the tissue 
m@auses electrical energy to be converted into heat thus 
raising the tissue temperature. A problem that is sometimes 
encountered is uncontrollable boiling of the tissue near the 
mome tip. A theoretical model which considers the combined 
emerects Of condliction, blood flow, metabolism, and rf heat- 
iMmeerts presented to show the relative effects of each of the 
terms on the temperature field. It was found that conduc- 
mio CfTrects and rr heating E€rrects aré Very important pa- 
maiierers in predicting lesion size, while the blood flow 
petect 15S Only marginally important. It is proposed that by 
Maimtaining the temperature of the probe at a low Constant 
fame, the maximum temperature of the lesion can be kept well 
below the boiling point, yet because currents and time can be 
Micheased, the volume of tissue destroyed ¢cquals or even Sur- 
messes that of the conventional probe. Finally, the radio- 
imeduency probes are compared to resistively heated proves. 
The lesions produced by these probes are shown to be signif- 


icantly smaller than the lesions produced by the rf probes. 
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1 eee Re DUGT TON 


Ns STATEMENT OF PROBLEM 

MIMOndehmbOntbhcat= Certain disorders Occurring im the 
central nervous system such as Parkinsonism, unbearable 
fom, and malignant tumors, surgeons at times find it neces- 
Sary to surgically destroy discrete regions Gi Mer yous cms. 
sue. Numerous methods have been used by clinicians to 
mecomplish tissue destruction. One such method, which has 
met with considerable SUCCESS, Cmpo 7s ce One a ClassiemnCec alec 
like probe which can be made to emit a low-power, continuous 
wave radio-frequency current. The probe itself serves as 


the emitting electrode and a second, indifferent electrode 


Mamet comas GOnmS remeee, sternal point om thie body aecic 
as a ground. As current is passed between the electrodes, 
tiem resistivity of the tissue causes electrical energy to 
moemeonverted to tissue internal energy thus giving rise to 
fmarencase 311 local tissue temperature. A lesion, or re- 
pero OL destroyed tissue, 1s formed when the local tissue 
Bemierature 1s raised above approximately 55°C [3]. 

The major drawback of the radio-frequency probe has 
Ween excessive heating of the tissue in the immediate vicin- 
mi OL the active electrode resulting in gas formation, car- 
bonization of the tissue and adherence of the coagulum to 
the electrode. Gas formation is many times accompanied by 


minor ‘explosions’ so that the method often reduces to a 


Bete. Drutal form of mechanical tissue destruction |35]- 


MZ 





Pimonech MLO Mimo cOnsmins predicting and creating 
well controlled lesions, it is necessary to model the heat 
MelolLer process OCCUrring in the region surrounding the 
merave electrode. Several attempts have previously been 
made at formulating such a model [3, 35]; however, solutions 
which consider the combined interaction of tissue conductive 
effects, blood flow effects, metabolic effects, and radio- 
frequency heating effects on the change of tissue tempera- 
ere with time have not yet been put forth. Such a model 
memoresented in this thesis for spherical and cylindrical 
Meowoes. the resulting equations are solved using a finite 
difference technique and presented in a form a can be 
readily used by surgeons. 

ee COUCSCH MLE HINGGMmOMN Une yCilmiine Ciee ag cm oo eee 
fieethe tissue, a constant temperature radio-frequency probe 
and a high resistance heating probe are also analyzed and 
amscussed. 

peaaditTona ll Ob)|CGtive 1S §hO presci Gears ah snot sence 
various techniques that are used to create lesions in t1s- 
me. tIhese techniques nay be broken into several broad 
Categories: mechanical, chemical, Morte radi1oO-acel ver 
X-ray, freezing and electrical methods of destruction. A 
description of these methods, complete with the advantages 


and disadvantages of each, follows. 


ies 





B. METHODS OF PRODUCING LESIONS 
ie echanical Destruction 

Perhaps the most conventional and familiar surgical 
technique is that in which the surgeon, through use of a 
scapel, manually makes an incision and removes or destroys 
unwanted tissue. This technique has been refined for use 
in the brain. A wire-like knife with a curved or angulated 
tip is concealed in a suitable length of hypodermic needle 
tubing of approximately 1% mm. diameter. The needle is in- 
e-otcd in the brain to the required depth, and the knife is 
pushed out and rotated 360°, thus producing a lesion in the 
mom Of a body of revolution. Glees and others [17] have 
used this technique and note that the lesions produced are 
foe Vocalim@ee but many times are large, irresulariy shaped 
peee lack uniformity. In addition, Carpenter and Whittier, 
ere port that "The lesions appeared to be dependent upon 
irecrruption of local vascular patterns and to a variable 
Peeeent tpon the blood elements liberated into the tissue." 
memo result of this dependence, the “rotating knife" method 
fae mot proven successful in producing consistent, uniform 
focaOns . 

More obvious disadvantages of this mechanical method 
mrerude profuse hemorrhaging, incidental tissue destruction, 
and infection. 

ee echem Leal Des krue clon 
‘localazed MeCrositse (lissic €eCStmuct ionic eprocd ea 


mwemically by injecting a few drops of zine chloride, chromic 
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Motion stilt lar Corrosave Chemical locally into the brain. 
iere Chemical tends to diffuse rapidly- and large, irregularly 
mmaped arces of tissue destruction are produced, Hyper- 
emeder lity, Convulsions, and unusual painware reported 

to follow lesions produced by this method. In addition, 
mmpoeted Solutions tend to “leak"™ out around the needle and 
@umse Cnlarged needle tracts and occasional cortical injury 
Bol] . 

BPeGalsce OL thew UnpredWerab lili ty Mor siicuSa ze ,meomaine 
m@emiaStOlogical characteristics of the lesions produced by 
corrosive chemicals, this method is of little value for use 
moeehne brain. 

eae Cenasonle Westruction 

Uist eso tc wames are mechanical vithieaitiens an 
solid, liquid or gaseous mediums lying above the range of 
Mmm nearing. The piezo-electric etfect 158 adapted to the 
generation of ultrasonic waves at frequencies above 500 
mmeoeyciles (kc), the range most often used in biological 
home. lissue 1s destroyed primarily by microscopic ee 
ical motion although heating 1s also present. Wood and 
moons |58| give an account of the effects of high fre- 
quency sound waves upon living matter. 

The ultrasonic waves are produced by a quartz 
meeate Vibrating in resonance with a high-frequency current 
Mmeecone through the plate. it is possible to focus these 
mmecceqate 2 Central focal point by using a concave quartz 


plate. One can then concentrate as much as 150 times as 


Ie 





nue@ieultrasoni1eG enerey at the tocal point as is found in a 
region near the vibrating plate [20]. As a result, maximum 
frssuce destruction occurs at the focal point while the in- 
tervening tissue remains relatively undamaged. 

ExMciiiicomes sp yea andmounpess | 20 | eshowsunan the 
ultrasonic technique causes an area of destruction on the 
scalp, followed by a region of normal, undamaged tissue, 
Piceernallysa well defined focal point of severe destruction. 
Winey also found that once a lesion is initiated, it appears 
m@emexeend rapidly and progress to non-irradiated parts of 
mae Drain. 

Presently (INS SUCCESS Ui WSS eb ne ultrasonic 
frenmique 1S hindered by the scalp injury. Methods are now 
meomo Styudaed to elamainate thas surface destruction. 

we noclodctive West ruction 

By implanting a small glass radon seed containing 
approximately 1 millicurie (mc) into a region of nervous 
mmosute if the brain, one can Create a sharp, circumscribed 
muearotf destroyed tissue of approximately 4 mm in diameter. 
Information concerning eidon dosage values for the produc- 
meme or appropriate neural lesions is included in a paper 
by Borison and Wang [4]. 

Stomp endmretersonmos | repOmtNtnac themes toms 
mmeduced by radon seeds are circular and discrete; however, 
they are variable in size and volume of tissue destroyed. 
Mey postulated that “The volumes of tissue destroyed were 


MomnopOmilonal te total radiation (malticuraes-hours ) > 
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pmeoblably Decause of Watterences in radiosensitivity of dif- 
imegent parts of the brain.” 

NGiweiougch data Were meaken to develop Conclusions. 
meee the experimental lesions produced were not lethal and 
did not result in physiologic disturbances. The lesions 
were also well localized and consistent in configuration 
and histologic characteristics. Perhaps the major drawback 
of this technique at present is a lack of radiosensitivity 
mmrormation for the various locations in the brain. A 
knowledge of this property is needed to effectively predict 
Mesion size. 

Toe ay WeStruct Lon 

By stereotaxically orienting multiple x-ray beams 
oa] , 6 15 8@heccrcticals yee ot omiaeias 
the intersection of the beams at a predetermined location 
meethe brain. This will result in tissue destruction at 
the focal point in much the same manner as the previously 
@mscussed ultrasonic technique. 

Carpenter ae Whittier [6] point out," Thesuscsor 
X-ray in stereotaxic technique appears promising, but such 
werhboOrate equipment in neuroanatomic laboratories is rare." 
In order to use this method, a detailed study would have to 
be under-taken so that stereotaxic orientation by mere loca- 
tion on the skull could be accomplished. 

'The advantage of this technique would be that pre- 
dictable and controlled lesions could be made in the brain 


PebcivoUt opening the skull. Presently, the therapeutic peams 


Jy 








Pe ailaple tor this method produce excesSive scatter, and 
beams small enough to be appropriate in an experimental 
minal are wot of suriicient intensity. There are x-ray 
beams that would produce the necessary intensity without ex- 
mectve Scatter, but their use for stereotaxic purposes is 
mre cessarily elegant. 
6. Cryosurgery 

Cryosurgery has proven to be quite usetul for var- 
fiemseapplications in the brain, especially for surgery com- 
plicated by bleeding or patient dabs COM Cuan | 

The cryoprobe is a hollow tube connected to a 
source of liquid nitrogen. Jhe probe is insulated throughout 
f-esept at its hemispherical tip. By controlling the flow- 
fmaeeme. Liguid nitre the tin can be cooled as loweas 
mio CGC, Causing the tissue in the immediate VEGi ni ty 7O eerie 
mepmco freeze and thus be destroyed. The cryosurgical tech- 
nique is bloodless and patient discomfort is minimal. In 
addition, the surgeon can produce only temporary numbing or 
Sumeaestroy a region of tissue. By observing a patient's 
Reaeti1ons to the numbing “tte Camene Sureecon, can ascertain 
fener Or not the probe 1s in the proper location before 
Hewerine the temperature to the lethal level. 

AedescLiptlOnmon tie = Les ton a Omite dace mene Ulaimoo ule 
fmeezes 15 given by Cooper and Trezck {10, 12]. A method is 
miso presented whereby the lesion Size LO: probe met mongoee 
trary dimensions can be predicted. The rate at which the 
temperature field develops is also considered, and a cryo- 


Surgical atlas is proposed. Both steady state and transient 


it 





Situations are analyzed in their work. Another heat trans- 
a modcleor the weryoprope 1S given by Rutkin and Barish 
m2 |. 

Cooper, Gionino,~and Terry (9! describe the effects 
which contribute to lesion formation by freezing: the slow- 
-ing down of the blood flow, congested blcod vessels, dehydra- 
Gren and the resulting concentration of electrolytes, 
me sctallization and the consequent rupture of cellular men- 
Deanes . 

Lack of hemorrhaging during the operation is a 
marked advantage of cryogenic surgery. The small amount of 
ee ehacings that JOS Ceeire as Cie 8/0) “elel= neeienseel pas- 


fee, Or the probe through the brain tissue. However, one 


mares probe size {2 mm. or greater). 
fee Uairect Current (dé) 

It is possible to produce generally well localized 
lesions in the brain by applying a direct current between 
Mremclectrodes: an active electrode inserted into the proper 
mreatlon in the brain, and an inactive electrode, or ground, 
Meeeced On some other position in the aoe ING Ae seiShS est Vale 
@= the brain tissue causes electrical energy to be converted 
memneat, thus causing the local tissue temperature to be 
raised. When the tissue temperature is raised above 55°C, 
the tissue is destroyed. 

UMisPtechniquevwean pe used with the active ellectrede 


Eemvying as either a@ cathode or anode. When the electrode 


Ins. 





Metenas a Gathodeé, the Lesions formed are large and irregu- 
lar due to the amount of gas formed at the tip [6]. In ad- 
mecton, Mullan and others |{22]) found the cathodal lesions 
momoe hemorrhagic and hard)to reproduce. Similarly the 
anodal lesions were also found to be irregular at higher 
-currents. According to Mullan's work, using the electrode 
as) an anode, surgeons can produce small and discrete lesions 
provided they slowly increase the current up to a maximum 
of 2 milliamps (ma) or prolong the duration of heating up 
to thirty minutes. 

In general, however, lesions produced by the direct 
@mrrent method vary in $size, shape and uniformity. Watkins 
weepecontizims the irregularity of the d-c lesions and adds 
cmd forse In omens t ieee . 
memeaesO States that “D-c devices Sutter from an ellectrode 
polarization and have an additional disadvantage that the 
make-transients and break-transients stimulate muscle and 
Memeous tissue." The polarization of the electrode is an 
Mmeemtrolled variable of electrolysis and, along with the 
amount of gas formed at the electrode tip, is a cause for 
miemiack Of control over the extent of tissue destruction. 

The Poet Ciciemammeet CUP Tem tac e mmnUelcite Sy etelNCmeiiean tr 
is an advantage in one respect: the surgeon can be assured 
mere the active electrede is in the proper location in the 
bain by observing a patient's reactions to the stimulating 
Peeece: Gt the direct current. 

The disadvantage of the nerve stimulation is more 


Obvious. Other nerves, as well as those directly in the 
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Mesion, could be stimulated during an actual operation. 
forcthermore, should the current find rts way to the heart, 
meevould cause the heart to quiver instcad of pump. ‘In ef- 
tee. the Suraeeon would have created an electric chair ef- 
eet . 

An additional disadvantage of the dc method is that 
feecee Current tends to have preterred paths in the brain 
(paths of least resistivity) resuiting in oval shaped 
Mesions being produced. 

One final disadvantage or the direct current method 
immmoted by Herrero [18}: “Animals with dc lesions usually 
ime extensive scar Eassue presen ee MOEN around sand silt iy male 
Memelin the tissue.'' This scar tissue implies that the blood 
ects ALOMUS With TNE NELVGUS CiSsSUuG Has been ccstroye 
feemene blood vessels remained, the body would more easily 
ceeemmodate the lesion. 

One of the advantages of the dc method over most 
fame methods discussed previously is that the small elec- 
trode can be painlessly inserted and causes minimal damage 
menemwsurrounding tissue while being inserted and retracted. 
im addition, hemorrhaging is not a significant problem with 
imme ce method. 

pee Rac Prequencys( re) Current 

Mme radio-ftrequency technique proaquces esionse1n 
men the same manner as the direct current method. An ac- 
mye, needle-like electrode is inserted into the proper io- 


Paaionein che brains and an inactive eleetrode 1S/plecci gon 
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the skull or some other remote portion of the body to serve 
as a ground. <A low-power, continuous wave, radio-frequency 
Siirrent is passed between the electrodes. Tissue destruc- 
tion is accomplished by resistive heating of the tissue it- 
pelt. tf the tissue temperature is raised above 55°C, the 
frecuc 1S destroyed. The factors governing the amount of 
heat produced are the resistivity of the tissue, the size 
oemene Clectrode, the amount of current used, and the length 
@meecime the current iS applied. An inverse relation between 
eirent and time gives rise to the amount of energy. If not 
Preven current 1S UWSed, mo E1Ssue destruction Occurs for an 
mmnice amount of time. J£ too much Current is ae Walle = 
lent boiling occurs immediately upon switching on the cur- 
meme . 

bPewmevseand Others elzZo | sconmdbeccd a SEUdsanconme la imme 
Meer destruction with the frequency of the aiternating 
Muemcines (hey observed that low frequency current resulted 
mmesporadic or no tissue destruction at all. If over 160 
Meeewere applied at four (4) Ki lohnertz (Chintz), “Sparks emanated 
Mmeolmerne electrode, As the frequency 1S increased above 45- 
Kinz and into the radio-frequency range, reac destruction 
mpecceds readily. Rowland, et. al. [Z8] offer an explana- 
mbometor the lack of coagulation at low frequencies. They 
pilasest that coagulation is dependent on the amount of 
charge per pulse. A unit of charge is a millicoulomb (mcoul.) 
emi 1S equivalent to a milliamp (ma) times a second. The 


moral energy applied by an alternating Current 15 determined 


1 





Peete number of mecouls per pulse times the frequency of 
pulses per second times the length of time current is ap- 
merccercealin that ume size ot the lesion is dependent only 
Sethe current (mecoul. per pulse times frequency) and the 
me. ltheretore, a low frequency cCtirrent applied in a given 
fMiemnitc Of time has a greater charge per pulse (mcoul. per 
mese) than does an equal amount of high frequency current 
meee Same amount of time. This greater charge per pulse 
could possibly be what causes the sparks to emanate from 
mites electrode. 

bwereterence, to the radre-¢requency technique, 


Puemiow | 5| explains gthat "Por surgical purposes it is nec- 


merary to have: (a) precise contro: of the size of the 
mieeieyorume O01 Gestruccion; (by; preéecisc control ci the 


pay oteat position of the lesion in a relatively inaccessible 
irmat2On; and (c) precise control of the functional position 
Geeeme Lesion." 

Wee “Bice ele. Gi ch eC oicenoe Pee yee foe jelniic: seeu@laforo 1. 2) — 
Meme y apparatus presented in the literature [1, 3, 13, 14, 
poems 52. 54, 355, 36, 37). Each of these set-ups contain 
certain advantages and disadvantages. 

ExXeessive ledtineror Ene ers cue ImmMoG rately: so lle. 
rounding the electrode results with the apparatus proposed 
byeescenow [5]. The temperature exceeds 100°C and the tissue 
mrens tO boll as a result. In addition, there is gas Lor- 
Marion, Carbonization of the tissues and adherence of coagu- 


iim to the electrode. As the hard coagulum becomes cncrustea 


Lo 





on the electrode, the total diameter of the electrode is in- 
creased and much damage to the surrounding brain tissue re- 
mes as the probe 1s withdrawn from the brain. The gas 
mormations are actually minor "explosions" so that the 
methoa reduces to a rather brutal mechanical method. Steam 
maepibes are tormed around the tip of the electrode and 
Mmeock the path of the current. An audible popping sound is 
mimrndi cation to the surgeon that boaling 1s taking place. 
In addition to the obvious disadvantages that occur with 
pemerine, it was noted that patients in which boiling occur- 
menor Len become “‘somnolent and confused" [35]. 

In an Seteuot bOmpleVem te) Olli cart Onno ce Umm apie 
mame den Berg and Van Manen [35] designed their "Coagrader" 
fomene nermits controlied and predictabie iesions to be rormed 
mememin Sate limits for physiological and clinical applica- 
moms. In addition to controlling and predicting the size 


Sere Lesion, the surgeon must be certain of the physical 





Hocation of the active electrode before irreversible damage 
Mmemrnitiated. As a méans Of assuring proper probe location, 
mereerequency alternating current is passed between the 
electrodes and the patient's response is analyzed. Sias and 
Seodman [31] point out that lowe requeney Oeouila & ive meu we Ile 
Mes the same stimulating effect as does direct current, 
whereas radio-frequency current has no stimulatory effect at 
asl . 


AnGther means ine Wale a .sureeon can sence. tio  prowe 


Meeation 15 discussed py Brodkey, ev. al. [Sj]. Tew 21scove) ca 
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that a temperature range of 40-49°C can reversibly block 
Pemeous TUNCLION- seni ter positioning the electrode stereo- 
meomaeeally and raisang the temperature of the tissue sur- 
femmes: the electrode to the reversible range, a patient's 
jmrpolse Will indicate whether or not the probe is directly 
Seman eect. Ihe range of temperature in which reversibility 
occurs is a narrow one, however. In addition, measurement 
Cemene temperature of the tissue surrounding the electrode 
meamm@ot a Simple task. It is easy enough to find the tem- 
Pemature of the probe tip by means of a thermal measuring 
@evi1ce, but the tip temperature is not necessarily the true 
mies uc eee re Surrounding it. It was shown by Dieckmann 
gmgeothers [13] that during the first few moments of heating, 
faeem@ete CL increase in temperature is greater aft a short 
distance away from the tip than at the tip itself due to 
Somauetive heat losses to the probe. Even externally ap- 
peed thermocouples have been reported to distort the tem- 
perature field of the tissue [5]. It has been shown, however, 
jfmecie radio-frequency technique has the fine control nec- 
Caedry to successiully reer ee Neene Sl byenveen| oy lee 

One ypeoblem that has been evident with the r-f 
Pemapinent concerns a defect in probe insulation. kLarge, ir- 
regularly shaped lesions have been known to form at points 
mone the electrode tract due to current "leaks." Szekely 
and others [34] developed a stylet probe to minimize this 
pleak ." 

caw th pele de probe. the Very silal ll csaze (Ome oc gag 


probe insures that the overlying brain tissue remains undamaged. 


“> 





SPecece probe so Intmomucecaus lowly, vascular structures can 


be pushed aside instead of punctured.. 


Closure rapidly takes 


place as the probe is withdrawn, and discomfort to the 


Beti1ent 1S minimal [39]. 


Merci Menem prebewis 1nSermeed inethe brain, and the 


Meeation of the probe has been checked, appropriate control 


Semene physical parameters of time, tip exposure, and tem- 


Menature allow well-circumscribed, 


be made [36]. 


mime Lesion produced. 


Geproducib lomslkesions te 


IWeweup diameter 2lso affects the size of 


ome tc: So recoded tienm are GlusSelUsscauby 


Pemey and Edsall {26}. There are 


tion: (1) coalescense around the 
meewly as the lesion progresses) ; 
NOGewthat , Unlike stire dc 


mores ollow preferred paths in the 


(EMMIS “SeRue one eos elile = 
electrode; (2) uniform 
oie B 
Geo ae O Flo, 
esieones eneet © Chest onmaec— 


Draliy DiLe. tothe redo. 


ieee all directions resulting in a well-circumscribed 


Spherical lesion. 


meme al core of dense dead tissue 


the lesroneecomteinse enroe  7ones © a 


surrounded by a zone of 


damaged tissue and finally, a narrow zone of granulations 


teen clearly mark the border of the lesion [2]. 


Magkewe cs. tare 


Roane through a lesion are often times spared. 


[ZiU] Societe inele eee lteyerel Wes Ge lls 


"This 1s 


Mewetially true when the temperatures used at the tip are be- 


itew sO°C." 
blood vessels. 


boiling temperature. 
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Heatine above 100°C results in rupturing tine 


hiomiacemoted earlier. thts Comme. DONG LO wuine 





Di outa ee radqto-ftrequency method allows con- 
meolled and predictable spherical HecTOns eo be produced 
weenin Safe laimits for physiological and clinical applica- 
tions. Hemorrhaging is not a problem and discomfort to the 
eeerent 1S minimal. There are two methods available with 
wimeh to check the probe location before irreversible des- 
fume? On Of the tissue begins. The system has precise con- 
meet over the volume of destruction and often times the 
blood vessels running through the lesion can be spared. 
Pemeaineg tends to be a problem with this technique, but 
proper control of various physical parameters can avoid 


boiling. 


ei 





ieee oNCAL ANALYSIS 


woe OUsmmlsctc ii the brain 1S quite com- 
meroewes) et Contains blood vessels, ventricles and is nonuni- 
Pm in Composition. Ihe grey matter in the brain seems to 
ifmawe a higher resistance to temperature than does the white 
Meeeet jis), Although the thermal properties (heat capacity, 
mrenmal conductivity, density and thermal diffusivity) are 
flewi tO Within an accuracy of 5% [J1], there are still many 
temmables in the brain whose exact values and temperature 
Sewitdeients are uncertain, 1.€., specific resistance, di- 
Emeernic Constant and circulation. 

Realizing these difficulties, a simplified, one dimen- 
Si0fial lticat transfer Wodei GL tue Yadsvo-LCLreyuency probe wiil 
now be developed. Consider the brain to be an infinite, 
memoeeneous, imperfect dielectric of resistivity p* and di- 
meeeeertc CONStant €. the active electrode 1s Considered to 
be a conductive sphere of radius r, embedded in the brain, 
while the indifferent electrode TS eu aesPpie nue dia e@ommane ton 
Aeenounding the brain (see Figure 1). 

The impedance of a spherical shell of material of width 
Ar can be represented by an incremental capacitance eh. Tih 


maeatlel with an incremental resistance dR [35] (see Figure 


Das: 





2 4nr* a eee 
dC = Ar > dk p Tie 7 
ee dR - dR 





et 





PROBE e4 WA 


MEDIUM INDIEFERENT 
ELECTRODE 


Figure 1. Model of Spherical Probe. 
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2 
dC = €° males 


Ar 
7 oe 
dR = 9 Aur 2 
_ dR 
We = 1+jwdkdC 
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The total impedance of the system can be obtained by in- 


moeomatame d’% {35}: 


At radio-frequencies, however, |jwp'e|<<1 for brain tissue 
ee) and the impedance reduces to Z = p'/4nr,. For ex- 


miele, tor brain tissue: 


[2 


Oe LSO= 200 cyanate [ay] 


2 


€ HOS al Sey) 


12. coulomb? 


e, = 8.85x107 PRS) y eee 





newton-m2 — ohm-cm 
e - 6.195x107°* SS 
Take 
o' = 200 ohm-cm 
fee n) 
o'e = 12x10°* sec. 
eniee , 
|jwo'e| = 12x10°+ w 
where 


Ges ini iz. 
Radio-frequencies range from 10 kHz to above the megahertz 
wea) £requency. Clearly, even in the megahertz range, 
Biop'el<<1, i.e, for LOMHz: 


10x10° eV 


Ijwote| = 12x107%<<1., 


WW 
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Pemaeresmuic, iethespadio-trequency range virtually all of 


Bi-meurrenespasses enmoucn tne resistive lege of the circuit 


shown in Figure 2. 


To solve for the temperature field emanating from the 


meano-Lrequency probe, a one-dimensional energy balance is 


fmepmeeed tO a Spherical ellément of tissue (see Figure 3): 


hector henepMemby BPCONCUICUCLTON comical Out or clement 
DPyeEecOnuct lon pied generadttonsdue toumetabolism + 
Heat input due to blood flow + Heat generation due CL) 
[mit eneatino -mimercaScurnieimntemmal eilerey On ole — 


Dine Te 


het ton (1) can be reduced by substituting the following 


- 


~~ 
i 


eS 





(i) CAO SO] SRC ELOnmel selec) So NS ie 
(ii) Elemental volume = AV = 4nr7Ar 

147 1 7 ] = = el Ol 2 oT 
(111) Heat in by conduction dan kA eS k 77 

; q din 
fey ) Heat out by conduction = den yy 

ele 
meh a hee aa 
Me or? 
(UES Oe Sees } 

(v) Heat generation due to metabolism = 


oe av) 
Smo¥ = Sy tet Ae 


iv 2) heeateaddatron due to blood tlow — 


mc, (T)-T)AV = mc, (T,-T) 4mr*Ar 


SW 








Sa (by “Conduction 





Heat generated 
aile. £O vie 
G@ligedsen) t 






Volumetric heat 
f 

generated due 

to metabolism 


(by conduction) 


Feeore S$. One-dimensional Energy Balance on Spherical 
Eb Lemeid ts 
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(vil) Heat generation due to rf current = 


coat 
2 _ I igeloluaine 
dk 4m r2 


(viii) Increase in internal energy = 


aT aT 


—ee = Z 
oc AV PC ae a eae Ay ee 


Note that in the model the effect of blood flow is treated 
Beeea peritect heat exchanger; that 1s, it 1s assumed that 

the blood acquires the tissue temperature as it exits a 
meee vOlLUMe Of tissue. This 1s consistent with a model 
fMmmcee set forth by Perl (24, 25]. Rosomoff and Holaday [27] 


give values for the volumetric blood flowrate My. 


If the metabolic heat generation rate is assumed to be 


wm ne 
te do ana A 


n mac tant Sf ran ha combined V7 LYTAw ta 
(ws WJ 4 WP wuss yg Re ae wer C4 4 A we wer Av ddd Dida NA ’ a m& Vb ke WY 


where 





This corresponds to redefining the reference temperature. 
ieeeas shown by Cooper and Trezek [10] that in the case of 
Page tissue, T, is on the order of 0.5°C higher than the 
hmeremic arterial blood temperature qT - FOr hie ease 2O bed 
radio-frequency probe, the blood flow convects a portion of 
the heat generated in the tissue due to radio-frequency ef - 
imeets away from the probe tip. 
EpolmsunstitueIne (1-Vijjd) aneo (1jeetie tollowime equa. 


e10n results: 
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_ 2 
Fr Ges Veh ts ~ 1) Aer ar 
(2) 
_ dT 2 2 
Be OG ar agin 2 2 ete hee 
10 
Mavadaing (2) by Ar and taking the limit as Ar>0 yields 
oid: 2 4 
in I je : De 0] 2 
pee geo, Me C, (T, Ey ae Oc ye 40 . (3) 
Mete that 
Erie oF 
fly RS ee 
meTICE , 


24 or 
a fara k aa F I2o! 


: aul , 
a2 Tee mop CE Me) Ai = oe Bye 


ibeek 1S assumed to be constant and (4) is divided by 47r°¢, 


imme) equation becomes 


k_ aa 


) — ear ee ee ere) ae (5) 


It is now convenient to introduce the following non- 


emmensional variables: 


2 
Me Cae 
bob 


R = — (Distance) 8 = es (Blood Flow Term) 
0 
Te on (Tame ) G = eee WS aa Heating Term) 
1 i ohiggeree kT, 
~ ie 6 
Se eer (Temperature) $¢ = &@ 


Baipstituting these variables into (5) yields 
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R= aa Se 7 B 6 Sue Re = oT (6) 


ihe boundary conditions ane 
ce eee prone mourttce w(r-1, Or R-1) no heat 15 eon- 
meered from the brain tissue back into the probe, or 


36 ALS 
aR |R=1 7 93 
(b) at large distances from the probe the temperature 


6f the tissue remains constant; that is, as ree, TT or 


0? 
imemnon-dimensional form, as R-7-, 670. 

icmuii tra becondstionm 15 thatvat tame t-0 1-1) or 6-0 
moeeeali locations in the brain. 

Mamvation (6) along with the boundary conditions and 
“initial condition can be DWtedn  taeites Clr therenee seam sana 
femeved)on Che digital computer. For the details of this 
M@eeractOn, the reader is referred to Appendix A. 


Pecimitar analysis €an be Carried out for a cylindrica: 


meaoee resulting in the following equation: 


Gr 2s ie : 
RoR °& * RF ~ Oy 7) 
where 
[2 ' 
G' = ees acai , on ae and B are the same as above. 


imesboundary conditions and the initial condition are the 
same as for the spherical probe. Equation (7) can also be 
Sollved by the finite difference scheme. The reader is referred 


Gemempendix B for the details of the cylindrical solution. 
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NOLC CimertOmupmen Ceomethies the rf heating term dies 
off rapidly as the radius increases. This implies that the 
iowmOrity Of the ri heating is generated close to the probe 
teeleD 

LuvasmMemaLonea in the Introdvetion, a problem that 
Peersts With conventional rf probes is tissue boiling near 
m@@emtip. fi the tissué temperature iS raised above approxi- 
mately 100°C, the rf technique reduces to no more than a 
eenae mechanical technique. It will be shown that even for 
iiea@erate Currents, unwanted boiling may occur. In order to 
Meewenmt Vaporization at the tip of the probe, it 1S proposed 
mmtenanecOonstant temperature probe be developed such that the 


maximum temperatures may be lowered and moved away from the 


« 


Tme SOLUcCi Ol CO this prObuiGin 15 Gecainced an vine 
mimrenmanner as tor the "insulated" sphere except that one 
bemmcdary Condition 1s changed, namely the temperature at the 
pee tip 1S now treated aS a constant. The results of this 
Memwod are discussed in the next section. 

nom me ie wecihmNqiuic used POmd@es ERO, micl VOUS 6C1 SS Uc mI Naene 
Memme employs a “high-resistance’ probe. In this method the 
M@imeent heats the probe tip rather than the tissue, resulting 
in heat diffusing through the tissue by conduction. This 
type of probe is discussed by Gildenberg [16]. One can 
solve for the temperature field by solving equation (6) with 
the rf heating term equal to zero and using the constant tem- 


Meewature boundary condition. Therefore, the equation to be 


solved is 


Sy) 








- BO = a. (8) 


This equation can also be easily solved using the 
finite difference scheme. ONCE - there is a closed form 
Solution to this type of problem stated in Carslaw and 
Miteger [7]. See Appendix C for details. This solution, 
Rieen Compared to the numerical solution, additionally serves 


mmonmcneck On the numerzecal scheme (Appendix C) . 
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Pils Wreoul loan pat SCUSS]ON 


Pipes EeviOUS atlLempc at moaclang the radio-frequency 
probe, Aronow [3] assumed that "Distribution of heat normally 
Mees place primarily by convection via the blood flow, 
Mmeener than simple conducticn."' He then postulated a model 
mre neglected blood flow and conduction entirely. It was 
mecumed that all the heat generated due to rf current went 
mieo ancreasing the internal energy of the tassue. The sim- 
fMeeetea model did not correspond well to experimental data. 
Papomow Indicated that neglecting conduction and blood flow 
muemuene MOSt probable reasons for this lack of agreement. 

Van Den Berg and Van Manen [35] have extended Aronow's 
model somewhat. Thev include the conduction term alongs with 
Mimreriro-trequency term but negleeg the blood tlow term azn 
waraolishing a steady state solution. They conclude that 
mienplood flow effect 1S negligible and that the model is a 
@eed approximation of the true system. Experiments with eg¢ 
white are in good agreement with their theoretical analysis. 

ms Was seen in the Physical Andlysis Secti0n, the pres- 
ent model includes the combined effects of rf heating, blood 
Mew. and conduction in analyzing the transient as well as 
Ge steady state solutions. Numerical results for both the 
emierical and cylindrical probe are given in Tables I-XxXV1. 
Values for 6/G are given at various non-dimensional locations 
(R) and.times (t) for various values of the blood flow param- 


meaeree (lables 124011). Similarly, values for 0/0, are7otvcn 
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miro rcONStant temperature Cylinder and sphere (Tables IX- 
Pav.) 

ie Pac ail Vener rcets (O. sconducetion, biood slow, time, 
proeeectirrent will now be discussed. In addition, the constant 
mimperature probe and the high-resistance probe will be com- 
Meee tO the insulated probe to show the relative advantages 


Go each. 


oe iii Peel Or CONvvCirONsOr THE TEMPERATURE FIELD 
Pucures + and 5S are graphs of the non-dimensional tem- 
memacuire versus the radial location in the tissue for the 
blood flow parameter, 8, equal to zero and non-dimensional 
meet Cduel to ten (corresponds to a real time of approxi- 
‘mately 2.5 seconds for a probe with a radius of 0.05 cm.). 
Mmaememirvyes are presented: the dashed line represents Aronow's 
memnerom which neglects conduction, and the solid line reépre- 
sents a solution which includes conductive effects. Note 
fmreretor a G of 1.0 (corresponds to applying a current of 
epemeuly greater than 2 ma for a 0.05 cm probe), 6 is 10 at 
Mmicmomnerical probe tip (R=1.0) if conduction is neglected 
(Aronow's assumption). This value is well in excess of the 
Boal ng temperature, ince POllane @oOcedirceat. ) 100 - Cee 
mimetial temperature of the tissue is approximately 37°C, 


therefore 


; _ boiling ~ be _ 100-37 | 
boiling Ie a 


its 4S Obvious trom Froures 4 and 5 thal conduction greae, 


reduces the temperature at the probe tip compared to what the 
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mipmeempecrature would be in the absence of conduction. 
Hence, a significant portion of the héat generated near the 


mmope £ip is rapidly conducted away from the tip. 


B. [tee brECleOr BLOOD FLOW ON THE TEMPERATURE FIELD 

Figure 6 is a graph of the temperature response in the 
Drain due to a radio-frequency current emitted from a spher- 
mel probe at a non-dimensional time of ten. Both conductive 
eeeects and blood flow effects are included in this graph. 
Peeail that all previous models assumed the blood flow ef- 
fect to be negligible. Note from the graph that the cooling 
Peace: Of blood flow decreases the temperature uniformly 
Mmimoeucnout the tissue. As the blood flow parameter, 8, is 
Mlereased, tre temperature 1S lowered further. Because ene 
blood passing through the lesion has a lower temperature than 
the heated tissue, heat is convected out of this region, thus 
reducing the paderature: TOP Cain taeewerc lane EC 1 Ste mit scaden 
meer ot this effect a typical value of 8 will be calculated. 
Pweeypical value for Mp Cy, Of  DLOOG- flow ile Gites wre 11) seee nel 
ey em’ /sec/°C [25] and a typical value of the thermal con- 
Guierivity of brain tissue is 0.0014 cal/cm/°C/sec PIs kere 
mmemm tissue, therefore, 8 1S appraximately equal to /7 Tr 
where r, is the probe fadius Expressed in Centimeters, “A 
typical radio-frequency probe might have a radius of 0.05 cm; 
mmeretore, 8 would be on the order of 0.02. As can be seen 
irom Figure 6, this value for blood flow marginally decreases 


mite eemperature. 
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Peewee Veisead eraphvotr the bleed flow erect assuming a 
aeemonical probe. Blood flow is seen to play a slightly 
meser role in Figure 7 than Figure 6. The higher values 
moines depicted on the graphs serve to illustrate the trend 
mrae the Dlood flow parameter has on the lesion size. Again, 
the blood flow parameter tends to decrease the temperature 
in the lesion resulting in a slightly smaller lesion. 

Figures 8 and 9 illustrate another effect of blood flow. 
Meeeenese figures the temperature at the probe tip (R=1.0) 
Mmeplotted against time for the spherical and cylindrical 
Probes. At low times the blood flow effect seems to be neg- 
MmiGiple, but aS time is increased, the effect of blood flow 
Beeomes Very noticeable. At t equal to 10 it can be seen 
mmouethe effect cf blood flow 
Maeeadius r, equal to 0.05 cm, t of 10 corresponds to approxi- 
mately 2.5 seconds. As the time is prolonged, the blood flow 
effect is increased even more. The larger values for 8 point 
out how increased blood flows cause steady state or equilib- 
fmm tO Occur at an earliér time. Again hinod eee is shown 
memaecrease the size of the lesion. 

Pamnce themprobe madiuswcemters the bleed tlow erneasee 
Samared term, one can see that as the probe radius 1s in- 
@eeased, the blood flow effects are sharply increased. For 
mpearecer probe the blood flow effect becomes a very important 
parameter. For the small probes and small times used in the 
Tf technique, however, the blood flow effects are only mar- 


Panally amportant. 








A) ab 
uo g JoJOWeIeg MOTA pootg sy} FO YI9FFI OvuL 


aqolg TeOtIputTAg e 


De Pesan aes 





46 








elegy Sa Ri ewiels Oat (1) See ym Si soit oc ay lee 
‘g ‘zoqyoWeIeg MOTI pootTg our FO (O'T=u) ATL 9qoig OY} 42 YD9FFT OUL “§ SANsTy 
, i/yn = 1 


haem 0°6 0°8 OZ 0°9 0°S O'b as 0°2 O°T 0° 0 


S00 0 


OTO°O 


SEO eanG, 





II 
Mm. 
© 
ee 
© 
a, 


0 


be] 





COGOLds (22 Pepa Ce LOL Pym OuT Tmsm@rte 4 











328 g SLOZPSOUPIRG MOT] POOT| 9y2 FO (O°T=y) dt oqold 043 2B 2OOFFA OUL °6 DANSTY 
3) ye <I) 

OO O¢ O'S OL 0°? OS OFF Ge. G4 O'/ 0°0 
— < «© gue ts? oes ey [on 
O° OT=9 

YY 
O°T=9 10 
70 
2'0 
"1 





SS 








ic. ieee? MeOrerNCREnvok OW Lite ON TEMPERATURE FIELD 

Diep ime ous secetion pointed out the fact that the ef- 
fect of blood flow _ miemeasedeas line Was, increased: 
mieure bE) amdicates the trend of the rf temperature field as 
the time is increased. One can see that at small times, the 
memo teliperatunre wmerease 1s large, but as times inerease, 
the rate of increase in temperature is lessened until steady 
mmec 15 reached. Steady state occurs at approximately T 
€qual to 150 (corresponds to about 36 seconds for a 0.05 cm 
maowe radits). 
Ie Peon Or rt HEATING TERM ION TEMPERATURE FIELD FOR 

iol 2D PROBE 

Pisircwl | 1hhustrabes Now tne pmebe tip Cenperaumre 1s 
Lewes qm len 
mee that a value of 6 of 1.7 corresponds approximately to 
MmewbOrlingy point. -One can see that a G of 5 is close to 
Gimesecritical value for 0 at the probe tip at time t=10. Jf 
MmitmresiStivity of brain tissuejis assumed to be about 150% 
miro, , the thermal conductivity cf the tissue is assumed 
to be 0.0014 Galen) -C/ scemun | the uli tially Cenperamire sO 
imemtrssue $/° C [3], and a probe of radius 0.05 cm, a G of 
mmeerresponds to a current of about 50 ma. In other words, 
memo cUrTrent Of 50 ma 1S applied for times longer than about 
pececonas, boiline will occur. o OLder EOumpreven Cubed dnice 


menuch lower current must be used, (about 30 ma). 
Be CONSTANT TEMPERATURE PROBE 
lt 1s proposea that boiling can be avoided altoBethpi if 


@ecomstant temperature probe is used rather than an insulated 
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Peobe ume touGe t2)isead graph dépicting the temperature pro- 
melee the brain for various values of the radio-frequency 
Meating term G for @ Eonstant temperature spherical probe 
Metco ©, B-0.1, t=10-0). As G is increased, the maximum 
temperature is moved away from the probe tip. The tempera- 
jmeeceprotile 1s more flattened than it was with the insulated 
mpeooe (see Figure ll). Note that the maximum temperature, 
aeweesor a G as high as 10, 1s well below the boiling point 
mee /), so there iS a great margin of safety with this 
Mien. Figure 135 gives a better indication of this safety 
m@etor. in this graph the insulated probe is compared to 
the constant temperature DEObemrer an ld Velu cl Uicwicm G sana 
feeeeiiie insulated probe is close to the boiling point for tT 


See whereas the constant temperature 


Ne « 


c nraltth hatau 
— o> tee ate od. 


tr ww a ft hs 


Seemporling point. The lesion radius depicted by the in- 


fmimeaeed probe for 7 of 10 is approximately 0.16 cm, i.e. 


Decrees Se Oran = oe 


9 a We eectonr ~ Eo eS Se eo 
besion IL; 47 


Peotlericure 13, 6=0.49E>>R=r/r,73.2. For r,=0.05 cm, r= 

wero cm or the lesion extends (.16-.05)cm cue. 1lecm £romeene 
meee For G=5.0, p=0.1, and t=10, the size of the lesion 
produced by the constant temperature probe is considerably 
Pireerer  iOweve:, the Constant temperature probe 1s well 
below the boiling point, whereas the insulated probe is dan- 


moeomems 8closec.  Invaddition, the current and time Can easily 
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Pemimereascds Without Coming Close to the boiling tempera- 
mince hon example, £rom Figure 12, it is evident that the 
constant ee ice wrote One 100 and G-10.0) pro- 
duces a lesion of approximately the same size as one de- 
meted an Fagure 11 for the insulated probe (8=0.1, 1T=10.0 
mew Ga-5.0), but the maximum temperature reached by the con- 
meme temperature probe is only Sie Mila eereater.s than 70 .C 
whereas the maximum temperature reached by the insulated 
mmowe 15 over 92°C. 

Another important factor Can Demo Ds ci vicde fl Onmuine sic 
mmaieese khecall that if the tissue temperature is kept 
below 80°C (61.16), My e tne i MenVvOUsSMErsSsuce in the Lesion 
bes) caestroyed while the blood vessels running through the 


a4. 1 ara liiac AL {- chaun 
at . ot &e Co ws 24 WY Ww 7 ate —_— at 4h re 8 4d. 


ER owemopamcan ou,,  Nesemm@met for 
Mememene COnStant temperature probe, only Ee iveom lc WV Olio (eG sik 
Mmemmeectroyed, whereas for G greater than or equal 5 for the 
insulated probe, the blood vessels as well as the nervous 
meee are destroyed. 

Mee nic Mapoint  COumMOeem=un reeanrdce tO the Constant wen 
Meteeure probe: not only can the current and time be manip- 
Meeecd tO determine the size and maximum temperature of the 
eeeton, but the temperature at which the PROCS NeW team 


mien oe Varied. By properly combining these three variables, 


optimum results can be accomplished. 


rs DG “SekbolSTANCE PROBE 
The final method to be discussed for producing lesions 


Mimene = Drailm 1S the high-resistance or pure Conduction probe. 





woe pOlmlted OWUEmea@riter, this technique consists of pass- 
mroecUrrene tnirough ome electrode and using the resistance 
Siecne electrode to Cause the probe to increase its own tem- 
MeGature, Fresulting in meat being conducted away from the 
probe into the surrounding tissue. This heat causes a tem- 
Memature rise in the tissue and eventually tissue destruc- 
tien results. Figure 14 shows the resulting temperature 
field for various values of t as compared to the insulated 
paeme temperature field. The lesions produced by the high- 
fmesmstance probe are significantly smaller than those pro- 
ame@ed by the rf method, but it should be pointed out that 
lemc xf ees ans have an added amount of Gicemeel nts semen aoe: Re 


milcweclecirical properties of the region are not that well 
mown, whereas the high-resistance le 
Gemedeat all on these é€lectrical properties. 

| In addition, by properly controlling the temperature 
emeeecie Clectrode tip, boiling is not a problem with the 


megi-resistance method. In summary, this method produces 


Sate, predictable and very small lesions. 


Gr. COMPARISON OF MODEL WITH EXPERIMENTAL DATA 

An attempt was made to Breck the proposed model with 
meabemmo cxperimental data {1, 3, 5, 6, 13, 28, 35, 36, 37]. 
However, this attempt proved unsuccessful due to the fact 
mrenot all the essential parameters were presented in the 
mepertved data. In order to correlate the analytical results 
With experimental data, it is necessary to know the size or 


Mi mOnODC tne Clurnent, applied, length of application, 
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temperatures recorded and corresponding positions in the 
Doro wr AemMeast One Of these parameters was either missing 
Smestated unclearly _ the Repoyted ee Pimenca lda tal. 
imnmordcnmmucmttluStrate this tact, a sample calculation 
will now be given from the data collected by Aronow [3]. 
ititesdata Consists of a graph of temperature rise versus time 
of application of current, and a table which gives the maxi- 
mum temperature recorded, the distance from the electrode 
at which this maximum temperature was recorded, the power 
used, and the electrode size. In OMe WET, iby Sige te elieeiedhsy 
Sieeea whether the electrode size refers to the probe radius 
or diameter. In Ach ELoire ONEwiSh tussUre Of tires Current 
Merde in the experiment. Aronow states that there is an am- 


= 
meter ana can mnt 
daiw & WwW 4aa re iw & v A ek 


eee neowc ue: , betricecoes: | t 
Mepulate the values of current or voltage used. 
It will be assumed that the electrode size refers to 


Mmeemueone Tadius therperore, ~,=0.05 cm. A temperature rise 


@f 45°C is reported at a 0.10 cm distance from the probe 


and a time of 6 seconds. Therefore, the following data is 
icequ@r ced : 

ie eo Cm 

To Os See OES 0S 9 10) Sl Sete 


ie >) See > Se 2 


t 6 sec. 
feewill be assumed that 


i, = SoC 


200 ohm-cm 


0! 
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= 0-001 "em’tsec 

Sage 
P cm ? 

7 Gail 
nec 

ss cal 
Pe ao JoMOEe cm°Csec ~ 0.0016 £2 
oc (1 gm ; of col yes Sec 
=i} aLon © 


Recall the following non-dimensional variables: 





2 = 
= ene oS = —_—— — = So es | Seay weer I ere ca 


An attempt will now be made to substitute the previously 
mentioned values into the non-dimensional variables and solve 


Horm the current. 








ety ee OSes ae 
ve i 0.05 cm ~ 
iF = 1th : 
§ = 5 oe ke eet ~ 1.22 
2 
Ti ( Ol ae es Baas cit} 
= wa Urls 
eal 
[0.0014 cmasec } 
2 
a (0.0016 si (« sec | | 
ee NE ne ee 84 
ze (a0 Semis 


Be iN owt gs ant J 


Lah “\fG 16n*r,?kT, 


i) ae, ee ee 

: OM ae kl — J , 
p 

See his Onn on od eandms: tS only 0.018) the effect of 


Pieeod flow will be neglected in this example. Jf a larger 


Sy 





meObee or a Longeyetime were Usicaewune Dlo@a 2low effect would 
meayea MOre important role. 
Enomelap len l whomean imsulated sphere at R=3 and T=3.84, 


fmenwcan be found that’ c/G 07078. 


_ 6 _ 1,22 . 

o = “j7e mye aoe! 
Now, 
~~ initial Yahi: 


(15.6) (16) (m7) (.05cm)?(.0014 eS" — ) (37°C) 





tb yolt ; eta ; RS Ceuk 


2) gon amps? VOLTtT amp watt sec 








1 = WV 67x10°° amp * 


] = Soe calibare 


This value for the current appears to be low (Aronow 
mmeeec that a typical value for current in 200 ma). The ap- 
Meeent discrepancy could be due to the fact that the probe 
Segenis in fact the diameter, or the distance from the elec- 
mgoeies 1S Considered to be the distance from the center of 
mmemelectrode rather than from the outer radius of the elec- 
meede., It the current had been reported, and the vagueness 
in the above terms eliminated, the data would seem to be 


more userul. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

SIpecoonsomancudoles to proamec comtrolled and predictable 
Semerical lesions in the brain with little or no hemorrhaging 
and minimal patient discomfort oe Pitema ad io 1 neq ema y 
Memenmmsaie. However, Care must be taken to properly control 
Meenwanount Of Current and length of time the current 1S ap- 
Memedeas CXCESSIVe Current or time results in tissue boiling 
miomunmeontrolled destruction. 

Mimoncr tO umnedlTeCtelesiOnm Size, one mMuUSt Understand 
mile cieectsS Of conduction, blood flow, metabolism, and radio- 
maeguiency heating on the tissue temperature. 

Conduction was found to play a very important role in 
Memenmiiniing lesion size. Ihe effect of conduction is to 
aeoersbute the heat away from the probe tip and carry the 
Memeeenco the brain. Ihe conductive effects very definitely 
mumet be neglected in attempting to predict lesion size. 

line biood =tlow ettfect was found to be marginally in- 
Meenant in determining the temperature field in the brain. 
Peemiding the etfect of blood flow tends to slightly de- 
mecase the temperature continuously throughout the field. 
mis, tie blood flow convects heat away from the region 
ground the probe. As a result, steady state or equilibrium 
merecached at an earlier time than if blood flow is neglected. 
The effect of blood flow becomes more predominant as time is 
gncreased. Therefore, if heat is being applied for any sub- 


Peiulrleanount Of time, the erfect of blood flow is relt 


onl 





Hones coerongly. the blood flow parameter, 8, is directly 
Eeeportiomal to the radius squared. As a result, blood flow 
plays a more caent any role in larger probes than for smaller 
probes. For the probes, currents and times presently used 

in the radio-frequency technique, the blood flow effects re- 
Main marginally important. One would overestimate the 

Meomon Size slightly by neglecting these effects. 

As the time of heating 1s increased, the lesion size 
MmaeneasesS gf a decreasing rate until steady state is reached. 
For a spherical probe of radius 0.05 cm emitting a current 
of about 20 ma, steady state is reached in about 40 seconds. 
Pia larger probe or a greater current were used, steady 
Meedee times would be larger. 


[oa — Se re ee revs CNN + ale bo aa a 9 ee 1 hi 
Wm me cS lie }|«686Cake UCCCm Ll Cin Ge OO eas ats Be ne 


hiew Shemp e mas 
4 


al 
Mmees in the lesion, especially near the probe tip as the 
femjerity of heating occurs in this region. As mentioned 
aimerer, Excessive current levels result in tissue boiling, 
whereas if current levels are too low, no tissue destruction 
occurs. 

Peete REC Demawinemacmunce probe tlp Gould be maintained 
mec lOW enough constant temperature by some cooling mecha- 
nism, both the current and the time can be Suet 1 Caitd youn 
@eeased, and a controlled well-circumscribed, predictable 
lesion would result without the Let Sit OCS wove mip emeibate. slave 
maximum temperature may be moved away from the probe tip and 
kept well below the boiling point yet high enough to destroy 


mervous tissue. In addition, the maximum temperature can be 
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Mepiemlowsemouch so tet the Dlood vessels running through a 
lesion are spared while the diseased nervous tissue is still 
mestroyed. 

iiestinaleproberanalyzed in the thesis is the high-re- 
fescance prope. This probe 1s also kept at a constant tem- 
Memacuyre, bUL the currént is used te heat the probe rather 
Miva the tissue. As a result, the conductive effect is the 
Sie mechanism by which heat 1s introduced into the surrounding 
meee. ihe Lesions produced by this method are signifi- 
ety esmaller than with the ri method. In practice, lesions 
Mmmeocieed USing the high resistance probe should be easier to 
paged 1 ct Ray, tidce created with the rf probe since only the 


fmierical properties of the tissue need be known. 


Pe 1S recommended that experimental work be carried out 
meer both radio-frequency probes and high resistance probes 
Memeneck the validity of the theoretical models. A thorough 
Mmieenbatire Search revealed sparse data reported; and these 
G@eaeae proved to be insufficient to substantiate the model. 
Not enough eaenet ere are given when data are presented. In 
Guger to effectively compare Meeuit<. Le ismemeccssary to. in 
eeude the following parameters when recording data: size of 
mem rOve, CUrrent applied, time of application, temperatures 
meeorded and Corresponding positions in the brain. 

It is highly recommended that the constant temperature 


radio-frequency probe be studied further. In particular, 


O5 





simple means of cooling the probe tip without greatly in- 
mreasine tne wsazcwOr the tip must be found. This probe 


Doula appear to hold great promise in brain surgery. 
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Pe ewer init DPRBERENCE APPROXIMATION POR A SPHERICAL 
PIO 


Statement of non-dimensional equation with boundary 


memertions and initial condition: 








5 oes: 
1 ok - pe +S, = 28 (A.1) 
(@) 0) 8g 
JOR} R=1 ; 
(b) As Rem, 8-0 
(c) At tT=0, 86=0 
where 
pet ee MEET o” oe cme 
i k 8 
a Ey 7 aio. 6 
a I . MOTT, P= 


There seems to be no closed-form solution available for 
Pwemecetvati1ol SO am CxpliGit, central difference, finite dit- 
Memetce technique was used to solve for the temperature field. 
pueeeeeeussj10n Of finite difference representation of various 
memms 15 fiven by Ozisik [23]. He examines the finite dif- 
Mmememce approximation of the terms by using Taylor Series 
expansions. 

Consider a spherical coordinate system with a region 
Mmterlaced with a R,t network such that the coordinates of 
a nodal point (i,n) is given as R=iAR (non-dimensional radial 


iM@ertion) and t=nAr (non-dimensional time). The non-dimen- 


™ 


. 
— + 
~~ 


Sional temperature at time t=nAv and radial location : 
: n 
@derimcd 2s O.. 
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Referencing Ozisik [23] the following transformations 


can be made: 


»(e a0) (A.2) 
r,t aR) \on_,) Non oan 
een FF ARe ‘( : wiew © (1+ Tek 28a i: 


B8 = BO. (A. 3) 
€ _  € 
Ra Gary: (A. 4) 
glug n 
6. - 6. 
Oo | il . 
ac a CEOS) 


Piescmeubins these finite difference approximations into 


eee) yields - 











Seuvinge for O. 


ntl _ At ( ; a RAS ee Oa 
ae BR {( 1) om ihe t)e ” , . AR? pat} OF * Gary 








ma 

@) 

ct 

a 

— 

I 

ee 
- 


ee : Bee. 2 n G AT 
oF = nf Q tery + (144 )¢ 1. + . 2M-BAR nf + TidR)*° 


tic ctibmiitemweritonia Serve as a GONStraint On the preb- 
mom, A Stability Condition must be introduced to limit the 
Haeimum Size Of time steps that can be chosen for a fixed AR. 


Mites Stability criteria for this problem will be stated without 
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proof: 1 - 2M - BAR*M must be > 0 (the interested reader 
/merererred to Ozisik [23]). 
If 1 - 2M - BAR*M = 0 


then 


i 


Me ee SNK Sere 


Mecall that M = At/AR*. Choosing AR automatically specifies 


Mee i.e. At = MAR*. Hence, 


ntl _ 7 n (1 1\,n GMAR?2 
; ba /?4-1 7 itt ete} * Garys: 9) 


Recall from (A.1) that the boundary conditions are 


(a) aye 
mel Ret - 


Ces RS OU, 


Heenes the finite difference scheme the boundary conditions 


become 


1p a8 


PittL ised 2 or 80, = Bays (A.10) 
where i* corresponds to R=1.0=i7AR, 
n n 
aaa Le = 0 or Oy.y = Nyy es 
where i = i*, i¥+l1, ... ,N-1, N. 
fre initial condition becomes 
6. = 0.0. | (A.12) 


ole 


OF 





Summary : 


tie berential equation 


. 
n+l _ Laie dean GMAR?* 
ae { ( + oa : GF )eaf * THaRy*° Yel) 


Boundary conditions 


(Ga) eet eee 
ie] ike] 
(b) gn Se 


jl 7 Neds 


ine ral condition 


( 
(c) me =.0-0. 


Commuter Program tL was used to solve this problem. 
Seeneardason’s Technique [8] was used to test the accuracy of 
meeesOlutzon. Three numerical solutions usins different 
Pmezedinerements of AR were obtained. Because a central dif- 
Hemenmce approximation has a discretization error of the ordcr 
@fyAR* [23], one can plot the approximation to the tempera- 
ture versus the increment squared (AR?) for each of the .three 
mmenementS used and extrapolate the plot to give the exact 
solution. Figure A.1 shows Richardson's technique as applied 
memeenie Spherical program. An increment AR = 0.1 proved to be 
Bmeenvenient interval for the program to work with and was 
few within the limit of accuracy desired by this approximate 
eoeution. 

An additional problem that was analyzed for this thesis 
MEME Nor case Ol 4 GONstanuntcmperdture probe. Recall tae 


+ 


ens probe serves to lower the maximum temperature incurre< 
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MmmeneephObe sSUyirace ana thereby serves to prevent the tis- 
Meron erOllin@ sine maititecrential equation and initial 
Sondition remain the same as stated in (A.13). Boundary 
Pendition (b) also remains unchanged, but boundary condition 
fyeis no longer valid. The probe iS now considered to have 
a constant temperature rather than being thermally insulated, 
Oe 


Oss = DONS eeynic. (A.14) 


Computer Program 2 was used to solve this problem. 
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eee x Bs COMPLETE SOLUTION FOR THE CYLINDRICAL PROBE 


Statement of nee 

Poetical probe son radius "r,s embedded in tissue 
foeen 1S anitidaliy at a uniform temperature of Tee eeaenat 
Mmitemtl1ssue there 1s volumetric heat generation due to metab- 
sesh Of strength Sa iimedathiron. there 15 blood £lowines 
momouch the medium with a flowrate My » Se Gia emnie ab Ch and 
memoerature ae emule ncmCNCuUl=tO-dlameter ratio Ox 
the probe is large enough so that edge effects may be ne- 
glected and assume the tissue is homogeneous and infinite in 
excent. <A radio-frequency eworieee | lS epee wets) seietayi telove 
Smebe, and the tissue 1s heated due to the resistivity p' of 
Feige (ones sic) ye Tala nll yehietsd) eubic 
mitcmerssuc and that none of the heat is conducted back 
momeouenh the probe, 1.¢€. the probe is thermally insulated. 
peesomeassume that at an infinite distance from the probe, no 


Marten is conducted. Solve for the probe temperature assuming 


Semeatant properties and one-dimensional radial heat flow. 


Notation: 
Symbol Description Tyjonkeruly Wickes 
OL Tissue thermal diffusivity Cay sce 
Ec Tissue specific heat Gade omy — 
Cy Blood specific heat ean omy © 
dR PiCrchoned stl oatemresi stance ohm 
J Cuasteemuie amp S$ 
k liver Malgmeonaleriva by 2Ole Tissue eal) Cisse 
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Symbol Deseriptrom yao eae Una ts 


my Volumetric blood flowrate gm/cm?/sec 
o' Resistivity of tissue ohm- cm 

0 Density of tissue gm/cm° 

r Radial location in brain cm 

oe Probe radius cm 

p Volumetric metabolism rate eal cn see 
i ane SEC 

a Tissue temperature sl Oe 

un Temperature of blood aC 

I Initial tissue temperature =16 


0 


Non-dimensional groups: 














femperature Distance Time 
o> 
0 
6 = 7 R = rt T = a 
0 To ; 
Blood flowrate re heating term 
Z 
a3 hy! Gt = 16 
k Ay seer 


feeece da One-dimensional energy balance on an element of tis- 
oemsce Figure B.1)}: Heat into element by conduction - 
Heat out of element by conduction + Volumetric heat generation 


wen co metabolism + Heat input due to blood flow + heat gen- 


eration due to rf current = Increase in internal energy of 
eeeoment . 
tea in by conduction = q._ = eee =: -k(rA@AzZ) ot 
ji arg) or or 
2 
din : dain 





meat OUtC by Conduction = Cee oe NB Ge: ST yee Arts 


(lay lor series) 


i 





Wolumetric heat generation due to metabolism = 
Se ero AOA) 
m m 
meat input due to blood flow = 
my Cy (Ty -T) AV = mpc, (rA6ArdAZ) (T, -T) 


fie ceneration due to rf heating = 


2qR = 29'AY 
peek eth ek 2 


miterease Im internal energy of element = 
oT 
OC sq (rACArAzZ) 
fomoinine the above 


2 
—_ aor 
or TASAZ 





+ mM, Cy C1, -T) radaraz “i 5, (rAsAraz) 











Chee) 
oC). y) 
fen. _i'p! - ae Sec ee 
ae saa mC, (Ty iain, 5,rAGAZ = pc ve PAGAN. 
Peeall gq. = -krA@AZ ore 
in or 
0(krA6AZ a) I20! 
— 520 | sakaaT my Cy, (CT, -T) rA6dz =P S FAG AZ 
- oe at 
= pc a5 MAGA Lx GoKSy 
Meevidineg through by rAgAZ 
| OT 
a( kr oT) 2 
1 dr ILS foes 7) = ONE 
ror * Sp * yenetaze * MpSplTy 7) = ec gE. BD 





out (PY GOncdvet1 on) 
Heat generation/ i 
Giicwo: Tie. ay / 
current a / 
‘ | Velumetric heat 
| Mime CenerTation due 
-_ IE) MS 1 ove Ie SiN 
i, , 
pod / Pleat mie dics Lo 
[ / ae piooda flow 
7 
rAo 
q; (by @omaic tion) 
meure B,J). One-dimensional Peto yee lance sOlm auGy | tice ct 
ip IEevenelte 


74 





fomraecOnseant k (B24) reduces to 


k a(x leon aT 
a eae? pe, er OS se CCB. S) 
feernts point a further constraint will be placed on the sys- 
meme Prior to the insertion of the probe it will be assumed 
mime the temperature of the medium is constant. In addition, 
it will be assumed that metabolism and blood flow have been 
active long before the insertion of the probe. 


If this be the case, then 


S 
m 


mp Cy CT, -T) 
or 

Wes ry + Sf My op: 
iMers implies that the medium temperature is warmer than the 
fluid by an amount Sf My cy prior to the insertion of the 
Meme. At Steady State the fluid thus removes the heat added 
ipmeteive VOlumetric heat source. 


Picci ne wroOuat Lon | B: 5) an thie form 


| ( Sa I20! aT 
eo.) ODS Ub oy. TS) * yenpzazze ~ © ae - B-9) 


sljeeests substituting T, for (Ty, ‘ Sf my Cy) When this is 


done, Equation (B.6) reads 


2 
ioe ee ere 


+ mo (1,1) * seneraze = °° 5 


oie 


Wet AG = 2% and AZ = I to get 


hs 





i are ® 7 OT 
ar + my Cy CT 9-7) + ena = OC Ce (B.8) 
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ale 


Mews tO Cast (8.8) into non-dimensional form let 





ae - 
Rel. gs ae ee ee hoe ei eg 2 !So! 
ie k Tag We 4n212kT, 


pUbStituting the above non-dimensional quantities into Equa- 


tion (B.8) gives 


a oi i) SS (B.9) 


Pama nome.) tS Ehemecauathton that must be solved in 
order to describe the temperature of the tissue as a func- 
meomeor tine when a radio-frequency current is emitted from 
meprope embedded in a medium with a volumetric heat source 
and a fiowlnug fluid. in order to soive the probiem, bound- 
Ppeeconaitions and an initial condition must be stated. The 
feunadary and initial conditions are as follows: 

(a) None of the heat generated due to rf current in 

Vii moore ms MeCONGCKeLCdual oO tie) pr OD cup Opiemi. 


enero 2 IL ane elise) ee 


ou 


o> or * 0 


ie i ori 
CE NOm- almenstonal form 
at R= 1, 90/dR = 0; 

(b) as ree, the medium remains unaffected 
aS LO lean 


Or, In non-dimensional form 


as Rem, 650; 
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US3) EE RE TOT eS Oe laces be 
or, in non-dimensional form 


at t = 0, 6 = Q. 
summary 


Woeererential equation 


ele 
Of R xe 
1 OR - G' _ 36 
[oe CP ee Sas total 
Boundary conditions 
30 . 
ae lrea = ° ets 
Ge (Bele) 


mitcial condition 


at t= 0, @ = O. 
Solution 


me velt. COl@ral Gditterence, 1ini te dirterenec, ap- 
peoximation was used to solve for the temperature field. 
Mensicer a cylindrical coordinate system interlaced with a 
R,t network such that the non-dimensional temperature at a 
time t=nt and radial location R=iAR is defined as Ee where 
Pema 1-1", wet)... ,N-1, N, and 1* is the location 
where R=1.0. | 

LcucmuiiPom COmO mci 2s lethc L£Ollowiine trans .Ornakrons 
can be made: 


a(x 22) 
y2o = i DRY ke {(i-pblo? +(1szt)o% 20 (8.19 


(ite) Bata es | 


lig, 





Emsertuting these finite 


wen 10) Canes 








il Jha \yat0l ( Isa 2 n G! 
AR? { (, on, . vest etaat AR? °a * (TAR 
les 
= ie Giles) 
solving for on gives 
(B.19) 
\ 
ght Aaa od ten " (Non % (. Zat : Pele Peg Ae 
i AR2 Sy gi AS alee 5S ae Gy aa) ees ae ne, 2 
a 2 
bet M At /AR (B. 20) 
ios de aso eal me 2 n G'At 
O = M {( Ao", + (Aor + {i ZENA per TiAR)? 
Meemaintain the stability criteria [23], 
1 - 2M - BAR*M = 0 (B.21) 
then 
a —— Ene 
FINI 
Recall M = At/AR. Choosing AR automatically specifies At, 
i.e. At = MAR?. Hence, 
(Baz) 
Oe = ak ro Pe Pee: G'M 
sie nf 7) °i-1 s(egyoraye | 


= BO. (B.15) 

=e B. 16 

= TibR)2 (B.16) 
ae 

- aot (B.17) 


Cech eMecwapprOxImations into 
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Recall the boundary conditions from (B.11) and (B.12) 


(a) 69 " 
ORIR=1— 


(b) as Re, 670. 


iemicecic Limite difference scheme, the boundary conditions 


become 
gh n 
1 *+1 aL ar Ik n n 
TAR 0 Or Oey Oe] (B.24) 
yea ~ ON-2 = 0 or 0. = 9” (B.25) 
ZAR N-1 N+1 ; 
Petoeewiy- i*, 7°+1,..., N =- 1, N, and i* corresponds to 
ieee, the initial condition becomes 
Q. = 0.0 B 
9 = Wer (B. 26) 
Summary 
Pree rcntial equation 
(cis 2 7) 


Boundary conditions 


(cas eee fait 
Spey 7 Speed 

Co) Ae aa 
Ne ON+] 


igtial condition 


Cy 
i 


By -specifying 8, the blood flow factor, and AR, the in- 


eecement., Mocan be found and At can be found. Computer Prosran 


iy 





3 does this and also solves for the temperature at various 
times and locations. Richardson's technique [8] was used to 
meek the accuracy of Pie wmactale rt recure B.2 shows Rich- 
meconm 5 tecimigque as at applied to the cylindrical program. 
Pi oevleom ele spherteal probe, the constant temperature 

meumdary Condition was analyzed with the cylindrical probe. 
Mremonly diztterence £0 the problem as stated in the previous 
Summary iS boundary condition (a). This boundary condition 


becomes 
Chae CONSE alle. Gb. 25) 


Computer Program 4 was used to solve this problem. 


SQ 





“WeIsOld TESTADULTA) EetOy SOU LUD opm a Os Perl Py uc: Gucins t- 


(av) 
On0'0 $00 _0£0:0 §to'9 Oro =sI0D =H SG D 


a a a Ns A 





AS 
Ce, See, 
$e ee x 4'0 
Mp pt 79 
80 
pe ee ee a a ee ee 
oa 


ol 


8] 





PoaweNnvi x Gs  SOLUTTONTFOR THE HIGH- RESISTANCE Sine CAL 
HEATING | PROBE 


DEeeie rete levedula tron 


oo a(Rt ER) 
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- (Gr) 

Boundary conditions 
i a On an cmeeOmMs tart ) (e272) 
AS Roe, 6-0. (G23) 
lmmtral condi t1i0n 
hit = (, 9 = 0. Cee) 
Epeapancdrme (C. 1) eeives 

2 30 3°6 _ 30 

Cer se 


Using a change of variables technique as described in Carslaw 


and Jaeger [7], let @ = v/R. 

















Then, 
2 (fi) 

00 7 LS eer aca : 
oR aR Zale ANS oe 
eG mee. letow  ov 1 ey 
ee |e Senn ol ses oo!) 

or 
ect ee oN OAV 
a9R2 RF R2 QR R QR2 

Substituting these new variables into (C.5) yields 
av, 2 av, av _ 2av,1 d*v _ gv. 1 av 
er R 3 R28] R dR2 "ROR GT 
ers vee Lay -_ 
RoR? 8 RO ROT aa 





Multiplying through by R 


EN 
a fees (C.9) 
meomecarslaw and Jaeger |/1, p. 134, the solution to this 


equation is 


_ 1 vie! OD for (R-1) _ | . YE (R-1) 





2¥T 
evfe JAR) 5 yee (C.10) 
on 
or 
= e- VB (R-1) erfc (R= 1) - VBT 
By 2 Te 
G@ei1) 
: (rR. 
+ eV (R-1) erfc —. + VBC 
27 
This equation was solved using Computer Program 5. The re- 


pets are shown on Figure 14. 


The steady state solution can be obtained directly: 


2 = a : : 
eer - Ov = Ompv=cje.% (RN , cie¥8 (R wd 


mepeying boundary condition (C.3) 


As Reo, +0 =>C, = 0 
_ 1 VB (R-1) 


eoRE= 2, 8 = 6, “PC, = 8, 


ce 





VB (R-1) 


2 
K 


Be 

O 
fmcomiparison Of the high-resistance program and the finite 
meelerence program for the constant temperature probe (tak- 
ing G = 0) serves as an additional check on the finite dif- 
ference program. Figure C.1 shows this comparison. As the 
meaph depicts, the two solutions fall exactly on one another, 
mus €1iVing additional assurance to the previous programs. 

When comparing the lesion size of the high-resistance 
technique with the lesion size of the insulated rf technique, 
Me is desirable to compare the size at steady state values. 
Bmmporder LO do this, the rf probe must be fixed so that it 


mioeune Same tip temperature as a comparable sized hign- 


recictanrea nroahe bas eto 


lyre 4) gee 
—_ ww af OY & SP Ae Se OY i eg al = ww _— od . ~_ 2 —— he 


tip temperature can be attained provided a proper value is 
chosen for the radio-frequency term G. Recall that for the 


rf probe the problem is as follows: 


Differential equation 


Boundary conditions 


= Boy = 
At R= 1, yR 0. 
As Rem, 6-0. 
lt € eae CONdd tom 


ei = Oe Oa a0, 
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Neglect blood flow and assume steady state 





_ 2» os 
TR 2 
230 _G 
Rv oR =~ RT SY 
98 _G | “1 
aR “Rt 
= eile oe, : 
But at R= 1, 22 -0=>C, = -G 
oe. © ~\oF 
aR RD 
me 
8 = - opr t Rt So 


Emieas Ro>, 620 =>C, = 0 





ap) 

iH 

i 

+ 
7] OD 


Pere ice ee hs Gobo ted Versus Ro Once 6 1s chosen: 
omewcan find a value for G that will result in the steady 
eaee Cemperature of the two probes being the same at the 


feeone tip. 
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